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Quantum-Classical Correspondence for Isolated Systems of Interacting Particles: 
Localization and Ergodicity in Energy Space * 

F.M. Izrailev 
Instituto de Fisica, Universidad Autonoma de Puebla, Apartado Postal J-48, 72570, Puebla, Mexico 



Generic properties of the strength function (local density 
of states (LDOS)) and chaotic eigenstates are analyzed for 
isolated systems of interacting particles. Both random matrix 
models and dynamical systems are considered in the unique 
approach. Specific attention is paid to the quantum-classical 
correspondence for the form of the LDOS and eigenstates, 
and to the localization in the energy shell. New effect of the 
non-ergodicity of individual eigenstates in a deep semiclassical 
limit is briefly discussed. 

PACS numbers: 05.45.-a, 31.25.-v, 31.50.+W, 32.30.-r 



I. INTRODUCTION 

In recent years, the growing attention was paid to the 
so-called Quantum Chaos, (see e.g. |l| and references 
therein). Nowadays, this term is used in very different 
situations and often leads to a kind of confusion. To 
clarify the subject, one needs to remind that the origin 
of this term relates to quantum systems which, first, have 
the well-defined classical limit, and second, in this limit 
the corresponding classical system is assumed to man- 
ifest strong chaotic properties. One should stress that 
systems under consideration are dynamical ones (or, the 
same, there are no random parameters in their descrip- 
tion). Therefore, this term was used in relation to the 
problem of the quantum-classical correspondence for dy- 
namically chaotic systems. 

Later, it was discovered that distinctive properties of 
Quantum Chaos (properties of quantum systems with 
the dynamical chaos in the classical limit) are generic 
for many other physical systems. For this reason, on re- 
cent conferences and workshops on Quantum Chaos there 
were many talks devoted to specific problems in atomic 
and nuclear physics, molecular and solid state physics, 
optics and acoustics etc. 

In order to classify the subject, we suggest to use global 
terms Quantum Complexity and Wave Chaos. The first 
term refers quantum systems with complex behavior (or 
complex properties of spectrum and eigenstates), both 
with and without the classical limit (see Fig.l). There- 



fore, we stress that properties of complex quantum sys- 
tems may have generic features, although the mechanism 
can be different. In the case of Quantum Chaos the 
mechanism is closely related to the deterministic clas- 
sical chaos, in contrast to quantum systems without the 
classical limit, where the mechanism of complexity is ei- 
ther due to a disorder (Disordered Chaos) or due to dy- 
namical reasons (quite complex interactions, although 
dynamical ones). As for the Wave Chaos, it refers to the 
chaotic properties of classical systems described by wave 
equations which are similar to the quantum Schrodinger 
equation. Therefore, many properties of such classical 
systems have much in common with those of Quantum 
Chaos and Quantum Disorder. There are many physi- 
cal situations in electrodynamics, optics, acoustics, etc., 
where chaotic properties of systems are well described by 
the methods developed in the study of Quantum Com- 
plexity. 
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FIG. 1. 

Below we discuss generic features of complex quantum 
systems, both dynamical and random ones, paying the 
main attention to the structure of strength functions and 
chaotic eigenstates. In what follows we consider the mod- 
els of isolated systems of interacting particles which can 
be represented as a sum of the unperturbed Hamiltonian 
Hq and the interaction V, 

H = H Q + V. (1) 

Here Hq describes finite number of non-interacting parti- 
cles or quasi-particles. The term V stands for a two-body 
interaction between the particles, and it is assumed to 
be responsible for chaotic properties of the total Hamil- 
tonian H . For disordered models the interaction V is ran- 
dom by assumption, however, the two-body nature of an 
interaction leads to important restrictions for chaos, see 
below. For dynamical systems the mechanism of chaos 
is not obvious, however, it is known to be related to 
classical chaos, if the classical limit exists. As for dy- 
namical quantum systems without classical limit, this 



case is less studied, however, the knowledge of proper- 
ties of Quantum Chaos and Disordered Chaos helps very 
much. We would like to stress that both the Hamiltonian 
Hq and interaction V can be integrable (for dynamical 
systems), however, the total Hamiltonian H may reveal 
strong chaotic properties. 

Such a separation of a total Hamiltonian (nf) into two 
parts is common in the description of complex atoms and 
nuclei and is known as the "mean field approximation". 
The core of this approximation is to choose such a ba- 
sis in which regular part of the total Hamiltonian is ab- 
sorbed in Hq thus resulting in new single-particle states 
(quasi-particles). In contrast, the residual interaction V 
accumulates all other terms which can not be embedded 
into Hq due to their very complicated structure. We note 
that the choice of the mean field is not well-defined pro- 
cedure, however, one can expect that if the most regular 
features of a system are described by Hq , main results are 
not very sensitive to a specific choice of the mean field. 



II. RANDOM MATRIX MODELS 

In the limit case of a very strong and " chaotic" interac- 
tion V 3> Hq, it is naturally to assume that system may 
be approximately represented by a completely random 
matrix of a given symmetry, this idea was in the origin 
of the Random Matrix Theory (see, for example, M and 
references therein). Giving correct predictions for sta- 
tistical fluctuations of energy levels of complex physical 
systems, such as heavy nuclei and highly excited many- 
electron atoms, full random matrices are too general and 
can not describe global properties of physical systems 
such as dependencies on the energy, strength of interac- 
tion, number of interacting particles etc. It is interesting 
to note that one of the first random matrix models stud- 
ied in 1955 by Wigner ( j|) has more complicated form in 
comparison with full random matrices. In analogy with 
(|l|), it consists of two parts, one of which is diagonal 
matrix with increasing entries tj and another is a band 
matrix Va with random elements inside the band, 



■"*.j ~~ e j"ij 



Vu 



(2) 



In original papers H the "unperturbed spectrum" was 
taken in the form of "picked fence", 6j = jD, where D is 
the spacing between two close energies and j is a running 
integer number, however, more generic case corresponds 
to random values 6j with the mean spacing D, reordered 
in an increasing way (<5y is the delta- function) . As for 
off-diagonal matrix elements Vij , they are assumed to be 
random and independent variables inside the band \i — 
j\ < b, with the zero mean and given variance, < Vij — 
> and < Vfj >— V Q 2 . Outside the band, the matrix 
elements equal to zero. Thus, the control parameters of 
this model are the ratio Vq/D of a typical matrix element 



to the mean spacing, and the band-width b. As one can 
see, the first term in (0) corresponds to the "mean field" 
Hq and the fact that the interaction has a finite energy 
range is directly taken into account. 

In spite of impressive results of the RMT in many dif- 
ferent applications (see e.g. 0), standard random matrix 
ensembles suffer from the lack of information about the 
nature of interparticle interaction. For this reason a new 
kind of random models was suggested in || (see also the 
review [pi); for a long time they were forgotten and re- 
cently they have again attracted much attention in the 
connection with the onset of chaos due to inter-particle 
interaction (see, for example, Q and references therein). 
Standard model of this kind describes n interacting parti- 
cles which can occupy m single-particle states (orbitals), 
for Fermi-particles it has the form 



H = J2 e s a\a s + ^Y1 Vs ^ 



Here matrix elements V s . 



S3S4 ^s l a S 2 U ' s 3 U 'S i 
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of the perturbation 



V^stand for a two-body interaction (indices Si, 82,83, S4 
indicate initial (si,S3) and final (.82,34) single-particle 
states coupled by this interaction), and e s is the en- 
ergy of a single-particle state \s). Creation-annihilation 
operators a\ . and a s . define the many-particle basis 



|fc) 



, a\ |0) of non- interacting particles, in which 



the unperturbed Hamiltonian Hq — ^2 e s a>to- s is diago- 
nal. It is convenient to reorder this basis according to an 
increase of the unperturbed energy Ek = ^2 e s n s with 
an increase of the index k = 1, . . . , N . In the Two- Body 
Random Interaction (TBRI) model all matrix elements 
V slS2S3Si are assumed to be random independent vari- 
ables. In realistic applications such as complex atoms 
and nuclei, matrix elements of the two-body interaction 
are calculated directly by using a proper mean field ap- 
proximation, see, for example, |7L|JJ. 

The TBRI model is very useful for the study of the 
role of a two-body interaction, as well as for establishing 
generic properties of Quantum Complexity, such as the 
localization, onset of chaos and thermalization in closed 
systems. To compare with Wigner Band Random Ma- 
trices (WBRM) of the type (g), one should stress the 
following. Having globally a band-like structure, the 
matrix Hij defined by (0) has many zero off-diagonal 
matrix elements inside the band, even in the extreme 
case when there are no additional integrals of motions 
and dynamical constraints, therefore, when all single- 
particle states are coupled by the (two-body) interac- 
tion. For m>n» the ratio of the number of non-zero 
matrix elements Vij to the total number of matrix ele- 
ments is very small, therefore, the matrix turns out to 
be very sparse. Moreover, non-zero off-diagonal matrix 
elements H^ have been found |nj] to be slightly corre- 
lated, in spite of the complete randomness of two-body 
matrix elements. The latter fact results in specific corre- 
lations between chaotic compound states and may lead 



to non-statistical behavior of physical systems with com- 
pletely random two-body interaction even in the region 
where conventional statistical description is assumed to 
be valid, see details in 10 1. 



The TBRI model (p) allows both to study the con- 
ditions for the onset of chaos and thermalization in 
closed systems, and to relate chaotic properties of spec- 
tra and eigenstates in many-particle basis to the prop- 
erties of single-particle operators. One of important re- 
sults obtained recently |ll| in the frame of this model is 
the Anderson-like transition which occurs in the Hilbcrt 
space determined by many-particle states of Hq. The 
critical value V cr for this transition is determined by the 
density of states pf — d'J 1 of those basis states which are 
directly coupled be a two-body interaction. The point 
that the interaction should be compared not with the to- 
tal density of states po of Hq but with pf for the first 
time was noted in |12j when considering specific model 
of a rotating nuclei (see also discussion in (b| ) . 

When the interaction is very weak, Vq <C df, exact 
eigenstates are delta-like functions in the unperturbed 
basis, with a very small admixture of other components 
which can be found by the standard perturbation the- 
ory. With an increase of the interaction, the number 
of principal components N pc increases and can be very 
large, N pc >• 1. However, if the interaction is still not 
too strong, n~ 2 ^/dfD <C V <C df , Q , the eigenstates 
are sparse, with extremely large fluctuations of compo- 
nents. In order to have ergodic eigenstates which fill some 
energy range (see below), one needs to have the pertur- 
bation large enough, Vq ^> df (for a large number of par- 
ticles this transition is sharp and, in fact, one needs the 
weaker condition, Vq > df ). 



III. LDOS: BREIT-WIGNER REGION 

In order to characterize global properties of ergodic 
eigenstates, it is convenient to introduce the F— function 
F- n ' which gives the envelope of eigenstates, 



p(n) = (n) 

3 — W 3 ' 



(n) _ 

■1 v ' 



c 



(") 



Cj(E^) (4) 



Here C are components of exact eigenstates \n) of the 
total Hamiltonian H expressed in terms of basis states 
|j) of H 0l \n) = V ■ Cj \j), and the bar stands for the 
average inside small windows centered at j or n. We use 
here the notations which refer low indices to the basis 
states, and upper indices to the exact (compound) eigen- 
states. Thus, the structure of eigenstates is given by the 
dependence wj 1 on j for fixed values of n . 

On the other hand, if we fix the index j and explore the 
dependence w " on n, one can find how the unperturbed 
state \j) is coupled to other basis states due to interac- 
tion. The latter quantity is very important since it gives 



the information about the spread of the energy, initially 
concentrated in a specific basis state | j) , over other basis 
states, when switching on the interaction. The envelope 
of this function u>: in the energy representation is re- 
lated the strength function or local spectral density of 
states (LDOS) which can be defined as follows, 



W{E {m \ j) = V \C ( ; i) \ 2 5{E - E {n ~>) 



(5) 



The sum is taken over a number of eigenstates \n) chosen 
from a small energy window centered at the energy E ^ m > 
. One can see that this function W(E,j) is originated 

(n) 

from the same matrix w which gives the structure of 
eigenstates, therefore, the shapes of strength functions 
and eigenstates are related to each other. Normalized 
to the mean energy level spacing, the strength function 
W(E, j) determines an effective number N pc of principal 
components of compound states \n) which are present in 
the basis state |j). 

For the first time the form of the LDOS for WBRM 
(g) has been analyzed in Ref. ||. It was found that the 
form of the strength function essentially depends on one 
parameter q — Po fc ° only. Wigner analytically proved g] 
that for relatively strong perturbation, Vq 3> D, in the 
limit q <C 1 the form of the LDOS is the Lorentzian, 



W B w (E) = — 1 ? 



2tt E 2 



lp2 : 

4 1 BW 



E = E-Dj (6) 



which nowadays is known as the Breit- Wigner (BW) de- 
pendence. Here the energy E refers to the center of the 
distribution. The spreading width Tbw (half- width of the 
distribution m)) is given by the famous expression, 



T B w = 27rp V" 



(7) 



In other limit q ^$> 1 the influence of the unperturbed 
part Hq can be neglected and the shape of the LDOS 
tends to the shape of density of states of band random 
matrix V, which is known to be the semicircle. 

Recently, Wigner's results have been extended in |15| | 
to matrices H with a general form of V, when the vari- 
ance of off-diagonal matrix elements decreases smoothly 
with the distance r = \i — j\ from the principal diagonal. 
In this case the effective band size b is defined by the 
second moment of the envelope function f(r) . Another 
important generalization of the WBRM studied in |l^] , is 
an additional sparsity of the matrix V, which may mimic 
realistic Hamiltonians. 

Random matrix models of the type (0) are very useful 
for the understanding some important properties of the 
LDOS. Let us, first, write the condition for the LDOS to 
be of the BW form lUl, 



D < T B w < A 



A F = bD 



(8) 



The left part of this relation refers to the non- 
perturbative character of the interaction, according to 
which many of unperturbed basis states are strongly cou- 
pled by an interaction. On the other hand, the interac- 
tion should be not very strong, namely, the width Tew 
determined by Eq.(f7|), has to be less than the width A# 
of the interaction in energy representation. The latter 
condition is generic for systems with finite range of the 
interaction V. One should also stress that, strictly speak- 
ing, the BW form is not correct for such systems since 
its second moment diverges. As was shown in 0, outside 
the energy range \E\ > Ae the LDOS in the model (g) 
decreases faster than pure exponent. 

For other models, such as the TBRI model where the 
band width b of the interaction is not well-defined, in- 
stead of Ae it is more convenient to use another impor- 
tant quantity, the variance a 2 , of the LDOS. The latter 
quantity can be rigorously expressed through off-diagonal 
matrix elements of the interaction, Cg = ^ Vh for i ^ j, 
therefore, a 2 = 2bV 2 for the model (g). As a result, we 
have Ae — koq/Tbw and Eq.(|8|) can be written as 



D <C T B w <C croV^ 



(9) 



Numerical data Jl|,||,|§ for the WBRM and TBRI 
models show that on the border Tew ~ 2&0 the form of 
the LDOS is quite close to the Gaussian, and this transi- 
tion from the BW dependence to the Gaussian-like turns 
out to be quite sharp. One should note that the form of 
the LDOS determines the dynamics of wave packets in 
the energy representation, see the data for the WBRM 
in Ref. M. 



IV. LDOS: TRANSITION TO THE GAUSSIAN 

For a long time it was assumed that in real physical 
systems the LDOS has the universal BW-dependence m) . 
However, when studying the structure of the LDOS and 
eigenfunctions of the Ce atom, it was observed clear de- 
viation from the BW-shape. Moreover , recent numerical 
investigation of nuclear shell models || has shown that 
the form of the LDOS is much closer to the Gaussian 
rather that to the BW. This fact is due to a quite strong 
interaction V ~ Hq , since the mean field Hq often in- 
cludes a large regular part of the interaction, thus leaving 
a "disordered part" (residual interaction) in V . 

As one can see, it is of great importance to find an 
analytical description of the LDOS in dependence of 
the strength of interaction. Although the extreme limit 
of a very strong interaction, q 3> 1, (or, the same, 
Wbw > Co) has been studied by Wigner in the WBRM 
model (||) , the semicircle form of the LDOS seems to 
be unphysical. Indeed, this form is originated from the 
semicircle dependence of the total density pv(E) , and 
the latter is known to be an artifact of the standard RMT. 



As was shown quite long ago (see the review pfl and refer- 
ence therein) , this result formally corresponds to random 
n— body interaction between Fermi or Bose particles with 
n = 1, ...oo. On the other hand, physical interaction is 
typically of a two-body nature, this fact is directly taken 
into account in TBRI models of the type (pi). 

It is clear that analytical treatment of the TBRI mod- 
els is much more difficult in comparison with full and 
band random matrices. For this reason, numerical data 
are very important since they may give a hint for rigor- 
ous results to be proved analytically (see, for example, 
0). The form of the LDOS in the TBRI model (g) 
in the extreme case when Hq can be neglected, is also 
defined by the density of states pv associated with the 
random two-body interaction V . As was shown in W, 
for large number of particles and orbitals, m 3> n 3> 1 , 
the density of states py has the gaussian form. This re- 
sult is expected from the point of view of combinatorics, 
however, the rigorous proof is non-trivial. 

Using general result of the WBRM model, see Eq. ([|) , 
one can expect that if the spreading width Tew > deter- 
mined by the density p/ = all of directly coupled basis 
states, larger than the "width" 2ctq of the LDOS, 



A 



1 BW 

~2oV 



> 1 , Tew = 2irp f V 2 , V 2 = (V%) , (10) 



the form of the LDOS approaches the Gaussian. As is 
indicated above, the variance ctq of the LDOS is defined 
by the off-diagonal matrix elements of Vi^j and can be 
found analytically (see details in |18 14 1), 



N 



(°S), = £^ 



v 2 



n(n — \)(m — n)(m — n + 3) 



j=i 



(11) 



It turns out that for Fermi-particles the variance of the 
LDOS is independent of the index j (therefore, of the 
total energy E^> ), in contrast to the models with Bose- 
particles, for this reason in what follows we omit the in- 
dex j for <tq . One can see that with an increase of 
the interaction, the half-width T/ lw of the LDOS changes 
from the quadratic dependence Tf lw ss Tew ~ Vq 1° the 
linear one, Thw ~ Vq. Sometime, this fact is missing in 
the analysis of the TBRI models. It is important to stress 
that in the case of A > 1 the parameter Tew has nothing 
to do with the half-width Thw of the LDOS, the latter is 
proportional to Co (see also discussion in M). 

Very recently, the form of the LDOS for the TBRI 
model (pi) was analytically found |l9| in general form. 
Since the density of states p(E) strongly depends on 
the energy, one should take it into account explicitly, 
W k {E) = F(E k ,E)p{E) . Here, F(E k ,E) is the 
F— function discussed above, written in the energy rep- 
resentation, E = E^ is the total energy (energy of an 
exact state \i) of H ), and E k is the unperturbed energy 



(energy of the basis state \k) of Ho . The method used 
in |ll| , is an extension of the approach developed in [p0| , 
which takes into account specific structure of the Hamil- 
tonian H . The dependence of Wk(E) was found to have 
the form, 



W k (E) = 



1 



T k (E) 



2n(E k -Ey + \rl(E) 



(12) 



with some function Tk(E). In the case of a relatively 
weak interaction, 1 < A < 1 , the function Tk(E) is 
almost constant on the scale of the energy width A# , 
therefore, the conventional BW-expression (H) is recov- 
ered for Wk(E) . In contrast, in the case of strong inter- 
action, A > 1 , the dependence of Tf.(E) on the energy 
E can not be neglected and the nominator in (O) is the 
leading one. As was shown in p9| , the function Tk(E) 
is defined by the density of states for one-body and two- 
body transitions, and for large number of particles and 
orbitals has the Gaussian form with the variance <7o . 
Detailed numerical study (l6| of the form of the LDOS 
in the region Tbw > 2cr have shown that the LDOS 
coincides with the Gaussian with a very high accuracy 

It is important to stress that the condition for the exis- 
tence of a smooth energy dependence T^ (E) is defined by 
the condition Vo > df of ergodicity for the components 
C k of the LDOS (or, the same, for exact eigenstates) . 
Below this transition, for Vq < df , the smooth solution 
of the equation for the function Tk(E) does not exist. In 
this sense, the result obtained in |l9| can be treated as 
the independent proof of the onset of chaos in the TBRI 
model when Vq > df . 



Having formal solution (12) for the LDOS, it is, how- 
ever, convenient to find a simple phenomenological ex- 
pression which depends on two control parameters only, 
the width Tbw defined by (|l0|), and the width Co , see 
Eq. (|ll|) . This was recently done in jl(|] where both the 
WBRM and TBRI models have been used to compare 
numerical data with the following expression, 



W(z;\,f3) = A- 



exp 



h<) 



z l + 



4 l BW 



(13) 



Here z — (E — E c )/ao is the normalized energy which 



refers to the center of the LDOS, with ctq given by (11) 



Taking the parameter [3 as the fitting parameter, two 
other parameters A and 7 are determined from the nor- 
malization conditions J W(z) dz — 1 and J z 2 W(z) dz = 
1 . As one can see, in the normalized energy units there 
is only one control parameter A . In the region A <C 1 the 
BW-dependence for W(z) is recovered, and for A ^> 1 
the Gaussian emerges, see details in |l6| . 

Numerical data show that this dependence gives quite 
good description of the LDOS for any values of A defined 
by Tbw and Co . It turns out that the transition from 
the BW-dependence to the Gaussian is quite sharp and 



takes place when A rs 1.0 (see also recent numerical data 
pl[). The above phenomenological expression (113) with 
the fitting function /3(A) is quite useful in the applica- 
tions. 

So far, we have discussed the form of the LDOS, pro- 
vided the eigenstates are chaotic in the sense that the 
fluctuations of their components (7- around the smooth 
envelope given by W(E) are gaussian-like. Similar anal- 
ysis can be performed for the form w{E) of exact com- 
pound eigenstates in the unperturbed energy represen- 
tation. As was noted above, these two functions, W(E) 
and w{E) are related to each other through the matrix 

wf ] , see Eq.(g). Analytical estimates @ for the TBRI 
model show that for not very strong interaction, when 
cr <C AE with AE as the total width of the energy spec- 
trum of Ho , the two functions are close to each other. 
This fact allows to use the results obtained for the LDOS, 
when considering typical structure of exact eigenstates. 



V. DISTRIBUTION OF OCCUPATION 
NUMBERS 

The knowledge of the F— function is very important 
for the relation of chaotic properties of the Hamiltonian 
H in the many-body representation, with the properties 
of single-particle operators. Indeed, exact eigenstates \i) 
of the total Hamiltonian H are given as 



|*>=£C«|fc>, 



!*> = < 



|0) (14) 



where C k is the k 



th component of the compound 

state \i) in the unperturbed basis. These coefficients Cj_ 
determine the so-called occupation numbers n s , 



,W 



E 

k 



c, 



(i) 



(k\n s \k) 



(15) 



where h s = a\a s stands for the occupation number op- 
erator. This quantity gives the probability that one of 
n particles occupies an orbital s specified by the one- 
particle state \s) for the fixed exact state \i) . According 
to this expression, this probability can be found by pro- 
jecting the state \i) onto the basis of unperturbed states, 
for which the relation between the positions of all par- 
ticles in the single-particle basis and the specific many- 
particle basis state is known by the construction of the 
latter. One can see that the probability n s — n s (E^') 
is the sum of probabilities over a number of basis states 
which construct the exact state. For Fermi-particles the 
occupation number n s = (k\ n s \k) is equal to or 1 de- 
pending on whether any of the particles in the basis state 
\k) occupies or not the single-particle state \s) . 

For chaotic eigenstates the n s — distribution is a fluctu- 
ating function of the total energy E — E^ 1 ' of a system, 



due to strong (gaussian) fluctuations of the components 

(i) 

C k . fn order to obtain a smooth dependence, one should 
make an average over a small energy window centered at 
E^ % > , which is in the spirit of the conventional statistical 
mechanics for systems in the contact with thermostat. 
In fact, such an average is a kind of microcanonical av- 
eraging since it is done for the fixed total energy E of 
a system. Therefore, one can define the n s — distribution 
through the F — function, see (0) and (Uq), 



,(£)= £) F(£ fc , £?<*>) (k\n s \k). 



(16) 



The n s — distribution plays essential role in the statisti- 
cal approach to finite systems of interaction particles, see 
details in [Mjlq]. It is clear that non-trivial part in the 
above expression for n s is the F— function F(Ek,E^) 
which absorbs statistical effects of the two-body interac- 
tion V. The important point is that in order to find the 
distribution of occupation numbers, one needs to know 
the envelope of exact eigenstates in the basis (or energy) 
representation, not the eigenstates themselves. There- 
fore, the knowledge of the F— function and the density 
of states po(E) of the unperturbed Hamiltonian Hq (it 
appears when passing from the summation to the inte- 
gration in Eq.(p6|)) gives us the possibility to relate global 
properties of eigenstates to the distribution of occupation 
numbers of single-particle states. 

The form of n s in the TBRI model (H) has been stud- 
ied in detail in connection with the problem of thermal- 
ization and of the onset of Fermi-Dirac (FD) distribu- 
tion. It should be stressed that the occurrence of the 
FD-distribution in closed systems of finite number of in- 
teracting particles is far from being trivial, especially, 
when the number n of particles is not very large. It has 
to be noted that in many applications n is the number 
of particles (quasi-particles) above the Fermi level, there- 
fore, it can be relatively small. For example, in the mean 
field description of the Ce atom there are n = 4 interact- 
ing electrons in the outer shell J7],D, an d n — 12 in the 
standard s — d shell model of heavy nuclei p . 

One of the most interesting results obtained in the 
frame of the TBRI model (0) is that the FD-distribution 
may occur even for 4 interacting Fermi-particles, pro- 
vided the (random) interaction V is strong enough. In 
this case one can obtain an analytical expression for the 
temperature and chemical potential which stand in the 
standard Fermi-Dirac form for the occupation number 
distribution n s . Another important result is that even 
if the form of the n s — distribution is far from the FD- 
distribution, it can be described analytically, the sim- 
plest case is when the form of the F— function is close 
to the Gaussian. The study of the n s ~ distribution al- 
lows to determine the onset of thermalization in terms of 
the ergodicity of components of compound states in the 
energy range defined by the F— function. In fact, the er- 
godicity of eigenstates results in standard (gaussian-likc) 



fluctuations, both for components of eigenstates around 
the envelope (F— function) and for the n s numbers, when 
changing slightly the total energy of a system (see details 
in 1 18 14|||). The latter property (existence of a smooth 
dependence of the n s — distribution with gaussian fluctu- 
ations around this dependence) is , in essence, the exis- 
tence of a statistical equilibrium for interacting particles 
in a closed system. 



VI. QUANTUM-CLASSICAL 
CORRESPONDENCE 

One of the important questions is the quantum- 
classical correspondence for the shape of eigenstates (SE) 
and the LDOS. As was pointed out in Ref. p2J, there is 
quite simple approach for finding classical shape of eigen- 
states and classical LDOS. This approach is generic and 
can be used for different physical models p3] p5i . 

Let us start with the classical SE. We consider the 
total Hamiltonian in the form (Q) where the "unper- 
turbed part" Hq can be represented as the sum of single- 
particles Hamiltonians H® describing the motion of n 
non-interacting particles, Hq — ^Zfc=i ^fc(Pfe> Qk)- The 
interaction V between the particles is assumed to re- 
sult in a chaotic behavior of the system H . Now let us 
fix the total energy E of the Hamiltonian H (t) and find 
(numerically) the trajectory pk(t) , qk(t) by computing 
Hamiltonian equations. Since the total Hamiltonian H 
is chaotic, any choice of initial conditions Pfe(O) , Q'fe(O) 
gives the same result if one computes the trajectory for a 
sufficiently large time. When time is running, let us col- 
lect the values of the unperturbed Hamiltonian Hq (t) at 
fixed times t = T, 2T, 3T, ... , and construct the distribu- 
tion of energies Eo(t) along chaotic trajectory of the total 
Hamiltonian H . In this way one can obtain the distribu- 
tion w(Eq; E — const) . Comparing with the quantum 
model, one can see that this function w{Eq\ E — const) is 
the classical analog of the average shape of eigenstates in 
the energy representation. Indeed , any exact eigenstate 
corresponds to a fixed total energy E = const and it is 
represented in the unperturbed basis of Hq . Thus, one 
can expect that for chaotic eigenstates in a deep semi- 
classical region the two above quantities, classical and 
quantum ones, correspond to each other. 

In the same way one can consider the complimentary 
situation. Let us fix the unperturbed energy Eq and com- 
pute the trajectory p\, ' (t) , q^, ' (t) which belongs to the 
unperturbed Hamiltonian Hq (t) . Similar to the previous 
case, let us put this unperturbed trajectory into the to- 
tal Hamiltonian H(t) and collect the values of the to- 
tal energy E(t) along the unperturbed trajectory for dis- 
crete values of time. It allows us to find the distribution 
W(E; Eq = const) which now should be compared with 
the LDOS in the corresponding quantum model. How- 
ever, in this case one should perform an average over 



many initial conditions J3fc(0) , <7fc(0) with the same energy 
Eg , since the unperturbed Hamiltonian Hq is separable. 
One can see that the above two classical distributions 
w(Eq;E = const) and W(E;Eq = const) determine the 
ergodic measure of two energy shells, the first one, when 
projecting the phase space surface of H onto H , and 
the second one, when projecting the surface of Hq onto 



H (see discussion in 22 1). 



VII. TWO INTERACTING SPINS 

The first example is the model of two interacting rotors 
(see [£3|,|4| and references therein), 



H 



M, 



(17) 



with angular momentum L and M . This model may be 
used to describe the interaction of quasi-spins in nuclear 
physics or pseudo-spins in solid state systems. Compar- 
ing with Eq.(0), here H = L z + M z and V = L X M X . The 
constants of motion are H = E, L 2 and M 2 , which are 
connected by the relation E 2 < E 2 nax = (L 2 + 1) (M 2 + 1) 
for LM > 1. It is worth to mention that in ( J1T| ) the dy- 
namical variables L, M are not canonical, however, keep- 
ing L 2 and M 2 as constants, one can present H in the 
canonical variables. 

It is known (see references in [ p3p4| ]) that in the clas- 
sical limit the phase space consists of both regular and 
chaotic regions, the size of which depends on the values 
of L 2 and M 2 . We consider the simplest case L — M 
for which trajectories are regular when \E\ is close to 
Emax = L 2 + 1 , and are chaotic at the center of the 
energy band, E ~ . 

Quantization of this model can be done in the stan- 
dard way, with angular momenta quantized according to 
the relations L 2 = M 2 = H 1(1 + 1) where I is an integer 
number. Therefore, for a given I the Hamiltonian ma- 
trix is finite and the semiclassical limit corresponds to 
the limit I — ► oo and h — > while keeping L 2 constant. 
According to the approach discussed above, the Hamilto- 
nian is represented in the two-particles basis \l z , m z ) and 
has the dimension N = 2l 2 + 21 + 1. Due to the symme- 
try degeneracy with respect to the exchange of particles, 
there are symmetric and antisymmetric states which are 
not coupled by the interaction (only symmetric states are 
studied in HIH). 

As was pointed out when discussing the TBRI matrix 
model (|3|) , it is important to reorder the unperturbed ba- 
sis in such a way that the energy of many-particles states 
of Hq increases. This kind of sorting corresponds to the 
classical model and allows to establish quantum-classical 
correspondence for the shape of LDOS and eigenstates. 
After this reordering, the unperturbed Hamiltonian Hq 
has the shell structure, with the same value on energy in- 
side one shell. Diagonal matrix elements are given by the 



eigenvalues —2lh, (—21 + 2)Ti, . . . , 2lh of the unperturbed 
Hamiltonian Hq, and diagonal elements of the perturba- 
tion V vanish due to a particular form of the interaction. 

The global structure of the matrix H m ^ n is quite spe- 
cific, however, it has some similarity with the TBRI ma- 
trices. Apart from the principal diagonal Ha determined 
by the eigenvalues of Hq , the next (to the principal 
one) diagonals H n>n ±i correspond to transitions inside 
each _ffo-shell while the " arcs" connecting the two corners 
represent transitions between neighboring shells having 
AHq = ±2%. Therefore, the matrix Hij has only four 
off-diagonal "curves" along which non-zero off-diagonal 
elements are located. All other matrix elements vanish 
thus resulting in extremely sparse matrix H^ . In spite of 
a clear band structure, the Hamiltonian matrix Hij can 
not be compared with Wigner Band Random Matrices 
of the type (||) since its band size is not much less than 
the total size of the matrix. Indeed, the band width is 
b = 21 + 1 in the middle of the matrix, therefore, the 
control parameter b 2 /N is of the order 1 and the finitc- 
ness of the matrix is crucial (see details in |2^j2J]). It 
is also important to stress that non-zero matrix elements 
can not be treated as random ones since the model ( |17J ) is 
dynamical. All non-zero off-diagonal matrix elements are 
positive and the mean and variance of the distribution of 
these matrix elements depend on the classical parameter 
L 2 = h 2 l(l + l) only. 

Numerical data for the LDOS in the center of the en- 
ergy band , E rs , have been obtained in the dynamical 
approach discussed in the previous Section. More pre- 
cisely, this distribution was numerically calculated tak- 
ing a sample of chaotic trajectories having the same fixed 
values of E and L 2 — M 2 . Following these trajectories, 
one can calculate Ho(t) = L z (t) + M z (t) taken at equal 
instants of time and find the distribution of Hq. In this 
case the classical distribution can be also evaluated ana- 
lytically since the classical unperturbed Hamiltonian Hq 
is integrable. 

The result is presented in Fig. 2 in the normal and 
semi- log scales. The quantum distribution Eq.(H) was 
obtained by averaging over I + 1 values of individual de- 

(k) 

pendencies w t for the largest shell with eo = . Here, 
the normalized energies are used, eo = Eo/E max and 
e = -B/i? m ax , with Eq and E as the unperturbed and ex- 
act energies. One can see a remarkable correspondence 
between the quantum LDOS and its classical counter- 
part. The only difference can be found in long tails of 
the distribution, which manifests a kind of quantum tun- 
neling beyond the energy region given by the classical 
Hamiltonian. Also, one can notice some discrepancy at 
the center of the energy band, where quantum effects are 
also expected to be relatively strong due to a finite value 
of the semiclassical parameter I . 

A very good correspondence has been also found for 
the shape of exact eigenstates SE, see Fig. 3. Numerical 



procedure of the averaging was taken in the way similar 
to that for the LDOS, with the average over I + 1 exact 
eigenstates taken from the center of the energy band. 
Again, one can see some difference in the tails and at 
the center of the distribution, which is caused by pure 
quantum effects. 
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FIG. 2. The LDOS distribution (full line) and its classical 
counterpart (dashed line) for the case L — 3.5, I = 39, after 

[24]. 



speak about a kind of the correspondence. This means 
that for the dynamical model of two interacting spins typ- 
ical form of the LDOS is neither the Breit-Wigner nor the 
Gaussian. One should note that in this model there are 
no free parameters (apart from the total energy) which 
can change the relative strength of the interaction V . 

When comparing the shape of eigenstates with the 
form of the LDOS, see Figs. 2-3, one can see a strong dif- 
ference in the energy range associated with the width of 
energy shells. This fact is due to the difference in the 
density of states po(E) and p(E) for the unperturbed 
and total Hamiltonians, see details in p3 24 1. It was, 
however, shown that with a proper rescaling the LDOS 
distribution coincides with the F— function with a high 
accuracy. 



VIII. ONE PARTICLE IN A RIPPLED CHANNEL 

Now let us consider the motion of a particle in the 
quasi- ID rippled billiard [p6|, see Fig. 4. 



It should be pointed out that at the edges of the energy 
spectrum where the classical motion is regular, individ- 
ual eigenstates differ very strongly one from the other, 
therefore, one needs take a relatively large number of 
eigenstates for the average, in order to obtain a smooth 
dependence on the energy. This means that although the 
quantum-classical correspondence can be achieved in the 
regular region, it has no sense since the energy interval 
to be taken for the average can be of the order the whole 
spectrum. 



y = d + acos(x) 



v y 




FIG. 3. a) Shape of eigenfunctions (the F— function) in 
the energy representation (full line) and its classical analog 
(dashed line) for the parameters of Fig. 2, after [24]. 

The interesting question arises in the view of the form 
of the LDOS, see Fig. 2. One can suggest that this form 
is close to the Breit-Wigner dependence (||) which gener- 
ically occurs in random matrix models. However, the 
detailed study |2!|] rejected this possibility since only in 
a very small energy region close to the center one can 




FIG. 4. Geometry of the rippled channel, after [26] . 

It was shown [£7],£8| that the dynamics of a classical 
particle in this billiard manifests generic properties of the 
transition to chaos with an increase of the amplitude a 
for a fixed d. A possible experimental realization of this 
model with finite length is a mesoscopic electron wave 
guide. In p7| classical transport properties of this bil- 
liard have been related to its dynamical properties, thus 
giving a transport signature of chaos. The analysis of 
the quantum motion in the infinitely long rippled chan- 
nel with the periodic upper boundary y(x + L) — y(x) 
is useful for the understanding of universal features of 
electronic band structures of real crystals, propagation 
in periodic waveguides, quantum wires and films, see ref- 
erences in [£7|,[28|. Detailed study of the energy band 
structure of the quantum version of this periodic rippled 
channel is presented in Ref. |2£| . 

The control parameter that determines the degree of 
chaos in the classical model is K — ^- , and strong chaos 

occurs when K ~ 1. In the quantum description one 
should solve the Schrodinger equation H = \{j) 2 x + f>y) 
with the boundary conditions ip(x, y)— at y = and 
y = d + (icosx. Instead, in Ref. |27| , |28| the approach 
has been used which is based on the transition to the 



new coordinates u = x, and v — 



for which 
at 



d+acos(x) 

the boundary conditions become trivial, ^(u, v) 

v = 0, d. On the other hand, in the new variables the 

specific interaction between two degrees of freedom arises 

and the Hamiltonian acquires a much more complicated 

form, 



H = ^g-WPaffgWPpg- 



1/4 , a,0 = u,v, (18) 



which is simply the kinetic energy expressed in the co- 
variant form. The momentum is now given by P a = 
— iK[d a + jd a ln(g)] — —ifig^^-^dag 1 ^, where g al5 is the 
metric tensor, and g — Det{g a fi) = [1 + ecosw] 2 , e = a/d, 
(see |2£| for details). As a result, the original model of 
one particle in the rippled billiard is formally transformed 
into the model of two interacting particles. In this way, 
the complexity of the boundary in the original model is 
incorporated into the "interaction potential" . 

The solution for the wave function can be represented 
by expanding the function $fc(u, v) in a double Fourier 
series, ^ a {u,v;k) = Em=iE^L-oo CmnWtt".*) > 
where a stands for the eigenstate of energy E a (k), and 
</>„„(«, v) — ir~ 1 / 2 g~ 1 / i sin(mirv)e' l( - k+n } u and the Bloch 
number was fixed k = 0.1 in the numerical study. The 
index m refers to the mode (or channel) number and 
n stands for the Brillouin zone number. Note that 
<f>mn are the eigenstates of the unperturbed momenta 
squared (and, therefore, of the unperturbed Hamilto- 
nian): P^ k mn = U 2 (rrm) 2 cj> k mn , P^tn = U 2 {k + nf f mn - 
One should remind that each pair (m, n) is associated 
with one number I which labels the unperturbed two- 
particle basis and it is essential that the unperturbed 
basis of Ho has to be ordered in increasing energy, 
E°{k) = £ {{n + k) 2 + (nnr) 2 ). 

In the classical limit one should perform the coordinate 
transformation in the same way as it was done for the 
quantum model. The classical Hamiltonian has the form 
H = ^g al3 P a Pp, with the same boundary conditions in 
(u, v) coordinates. Expanding this expression, one can 
obtain, 



1 

H =2 



P z 



l + < 



PuPv 
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(1 + eO 2 °. 



(19) 



where £ = cosu. £ u = d£/du, and e = a/d. It should 
be noted that the transformation from the old to new 
variables is canonical. In accordance with the quantum 
model, the separation of the total Hamiltonian into two 
parts, H = H° + V , has to be done as follows, 



H ° = 1 [Pi + P 2 ] ; V=H-H° 



(20) 



Let us now compare the quantum LDOS and the 
shape of exact eigenstates with their classical counter- 
parts which are numerically obtained [g6[ in the described 
above approach, see Fig. 5. One can see, that, in general, 



there is quite good global correspondence between quan- 
tum and classical data (the role of quantum fluctuations 
around the envelopes will be discuss below). The inter- 
esting point is that the forms of the LDOS and the SE are 
highly non-trivial, if to compare with the Breit-Wigner or 
Gaussian distribution. Remarkably, the averaging proce- 
dure reveals a three-peak structure for both classical and 
quantum quantities. 
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FIG. 5. (above) Shape of eigenfuctions (the F— function) 
in the energy representation (broken line) and its classical 
counterpart (smooth curve); (below) the LDOS and its clas- 
sical counterpart. For the F— function the average is taken 
over exact eigenstates from the interval 1870 < a < 1950 and 
for the LDOS the average is over basis states from the interval 
1900 <j< 1990, after [26]. 

The origin of these peaks can be understood by a de- 
tailed analysis of the classical trajectories Pq| . Specifi- 
cally, it was found that even for the chaotic motion the 
particle can be located for a long time in the neighbor- 
hood of stable and unstable periodic orbits of period 1. 
One can find that the unstable (stable) periodic orbit is 
defined by P u = and x = 0; 7r respectively. The right 
peak corresponds to the motion perpendicular to the x 
axis at x = 7r (unstable orbit). Similarly, the left one 
results from the stable orbit at X = . In contrast, the 
central peak is originated from the trajectories which are 
nearly parallel to the x-axis. A similar analysis explains 
the origin of the structure of the classical LDOS . 

As for the comparison between the form of the LDOS 
and that of the SE, it is easy to detect a difference, es- 
pecially, in the central region. The reason for this is still 
not clear, in Ref. p6| it was suggested that it is due to 
the difference between time averages and phase space av- 
erages for finite times. 



IX. LOCALIZATION AND NON-ERGODICITY 

IN THE ENERGY SHELL 

Now let us discuss the problem of localization in the en- 
ergy shell. As was claimed above, with an increase of the 
interaction in the TBRI model, the Anderson-like transi- 
tion is expected [pTj in the Hilbcrt space of many-particle 
states. It occurs when the interaction strength Vq exceeds 
the mean level spacing between basis states which are di- 
rectly coupled by the two-body interaction, Vo > df . 
In terms of the structure of compound states, this means 
that above this threshold the components C ergodically 
fill an available energy shell. More specifically, the fluc- 
tuations of these components fluctuate around the mean 
given by the F— function, in accordance with the conven- 
tional statistical mechanics, typically, according to the 
Gaussian distribution. One should stress that for closed 
systems the width of the energy shell is always finite due 
to a finite range of interaction, and it can be much less 
than the total size of the energy spectrum. However, the 
number N pc of principal components can be very large. 
In this sense, one can rigorously define the onset of chaos 
(or maximal complexity in eigenstates), as the ergodic 
filling of the energy shell with the Gaussian fluctuations 
around the envelope given by the F— function. 

If the system has the classical limit, the F— function 
has its classical counterpart which can be easily found 
from Newton equations of motion. More difficult sit- 
uation is when there is no classical limit, in this case, 
although the F— function exists, it is of the problem 
to find it without direct diagonalization (and the aver- 
age over a number of eigenstates) of a (typically huge) 
Hamiltonian matrix. Using two examples of simple quan- 
tum dynamical models with the classical chaos, we have 
shown here that in a deep semiclassical limit there is 
a good correspondence for the form of both the LDOS 
and eigenstates. Strange enough, this quantum-classical 
correspondence is studied scarcely, there are mainly nu- 
merical results discussed above (see also, p5|]). As the 
data for the rippled channel show, a kind of semiclassi- 
cal periodic orbit theory can be developed both for the 
LDOS and SE. 

The knowledge of the SE is extremely important for 
the study of any kind of localization in the energy shell. 
The non-trivial role of the F— function is that it gives a 
reference to which individual eigenstates should be com- 
pared in order to see the localization. Let us first consider 
the TBRI model where there is no influence of classical 
periodic orbits which result in a specific sort of localiza- 
tion (scars, if to speak about the quantum-classical cor- 
respondence in the configuration space). For such kind 
of models, above the threshold Vq « df the gaussian fluc- 

(i) 

tuations for the coefficients C are expected. This was 
checked for excited states of the Ce atom which is known 
to be quite chaotic dynamical system PS|J7||. Numeri- 



cal data have shown, that, indeed, the fluctuations turns 
out to be very close to the Gaussian. Note, that for the 
study of these fluctuations, first, one needs to normal- 

(i) 

ize each value C- to its mean-square-root value which 
itself depends on the index i (or on the energy E"\ if 
to analyze the fluctuations in the energy space). There- 
fore, one should find this mean value (square root of the 
F— function) by a kind of the average. 

Below the threshold Vq ~ df the fluctuations of the 

(i) 

amplitudes C are extremely large since the components 
which correspond to those states, which are not coupled 
directly by a two-body interaction, are very small, thus 
resulting in a large peak in the distribution of P(C) 
around zero. As one can see, the quantum localization 
results in the sparsity of the F— function and the LDOS, 
in contrast to the standard Anderson localization in one 
and quasi- ID disordered models, where eigenstates are 
"dense" in the configuration space and occupy a finite 
number of sites (analog of basis states in our models). 
Therefore, the fingerprint of the Anderson localization is 
an exponential decay of the amplitudes Ck far from the 
"centers" of eigenstates, unlike the localization in the en- 
ergy shell for closed systems, where the cut-off of tails al- 
ways exists due to finitcness of the energy shell. On the 
other hand, for a very strong localization, when eigen- 
states consist of few basis states (for Vo <C df ), the 
situation is quite similar to that for the Anderson local- 
ization (this typically happens for low energy states due 
to a relatively large values of df ) . 

A much more delicate situation occurs for dynamical 
quantum systems with the classical limit. In this case in 
addition to the quantum localization, even in a strongly 
chaotic classical region, the influence of tiny structure 
of the classical phase space (for example, short unstable 
periodic orbits) may be quite important. In this sense, 
one can speak about the localization which is caused by 
both classical properties (periodic orbits) and quantum 
ones (scarring effect). For this reason, the mesoscopic 
region where both classical and quantum effects can be 
equally strong, is of the most difficult for the analysis. 

Coming back to the numerical data for the rippled 
channel Eq] , a noticeable fraction of extremely localized 
eigenstates has been detected in the region of a strong 
classical chaos where typically eigenstates are ergodic. 
It was found that these eigenstates are originated from 
those basis states which are located on the edges of the 
shells of Hq . It turns out that these basis states are very 
slightly coupled with the others by the interaction V . 
These localized states correspond to almost free motion 
of a particle along the channel. The presence of very lo- 
calized states has been also observed in the model (|17| ) 
of two interacting spin [t23|]24f | . In this case, their origin 
is not so clear and deserves an additional study. 

Finally, we would like to raise a new problem related 
to the quantum-classical correspondence for chaotic clas- 



sical systems. It refers to the ergodicity of individual 
eigenstates in a deep semiclassical region. According to 
the Shnirelman theorem, it is believed that the closer this 
limit, the less influence of quantum effects. To see this 
fact, a kind of the average over a small energy range is 
typically assumed. On the other hand, in a deep semi- 
classical limit the eigenstates have a very large number 
of components. This allows one to perform the statis- 
tical analysis inside individual eigenstates, without any 
energy or ensemble average. This is in a spirit of dynam- 
ical chaos for which statistical properties can be revealed 
along a single chaotic trajectory. Specifically, we can 
ask: how fast these fluctuations approach those ones pre- 
scribed by the standard statistical mechanics? Numerical 
data |2(| manifest that individual eigenstates can not be 
rigorously treated as random ones since non-statistical 
deviations seem to persist for any finite values of Ti. In 
this sense, one can speak about non-ergodic character 
of individual eigenstates, in spite of the existence of the 
classical limit itself. Definitely, this question deserves a 
detailed study, both numerically and analytically. 
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